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Abstract 

We find that quantum spin Hall (QSH) state can be obtained on a square-like or rectangular 
lattice, which is generalized from two-dimensional (2D) transition metal dichalcogenide (TMD) 
haeckelites. Band inversion is shown to be controled by hopping parameters and resnlts in Dirac 
cones with opposite or same vorticity when spin-orbit coupling (SOC) is not considered. Effec¬ 
tive k-p model has been constructed to show the merging or annihilation of these Dirac cones, 
supplemented with the intuitive pseudospin texture. Similar to graphene based honeycomb lattice 
system, the QSH insulator is driven by SOC, which opens band gap at the Dirac cones. We employ 
the center evolution of hybrid Wannier function from Wilson-loop method, as well as the direct 
integral of Berry curvature, to identify the Z 2 number. We hope our detailed analysis will stimulate 
further efforts in searching for QSH insulators in square or rectangular lattice, in addition to the 
graphene based honeycomb lattice. 
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I. INTRODUCTION 


Topological insulator (TI) has attracted significant attentions from researchers in re¬ 
cent years.^i®' Its bulk band structure is characterized by a global topological invariant Z 2 
number,!^ which is protected by time-reversal (TR) symmetry. The bulk-boundary corre¬ 
spondence ensures that on the boundary of a TI there exists topologically protected metallic 
states which cross the bulk band gap by connecting the valence and conduction bands. Es¬ 
pecially for two dimensional (2D) TI, the conducting one dimensional (ID) edge states have 
opposite velocity for opposite spin channel, which leads to quantum spin Hall (QSH) effect. 
Therefore, 2D TI is also known as QSH insulator. Since the backscattering of electrons 
in this ID edge state is prohibited as long as TR symmetry is conserved, dissipationless 
quantized electron conductivity is expected and has great potential application as an ideal 
conducting wire in device. 

However, up to now, only two materials have been confirmed by experimental observation 
of quantized conductivity inside of bulk band gap, namely, the quantum well structure 
composed of HgTe/CdT^^^and InAs/GaSb.^ There have been intensive efforts in searching 
for 2D QSH insulator. One is based on the honeycomb lattice, which is stimulated by Kane- 
Mele model for graphene.The low-buckled silicen^I^, chemically decorated single layer 
honeycomb lattice of Sn,^G^I^and Hi or StP^are all in this category. The second category is 
represented by ZrTes/HfTe^and Bi 4 BrJ^. They are single layer exfoliated from their three- 
dimensional (3D) counterparts, which are weakly bonded layered materials. The proposal 
of IT’ structure of single layer transition-metal dichalcogenide (TMD)ff^ can be ascribed to 
this type also. The third category are based on square or rectangular lattice. Such as recent 
proposals of 2D buckled square lattice BiF® and rectangular TMD haeckelites.^^®^ 

In this work, we will discuss the possible topologically nontrivial quantum states based 
on the model built for 2D monolayer transition metal dichalcogenide (TMD) haeckelites,^®^ 
among which MX 2 (M=Mo, W and X=S, Se, Te ) have been predicted to be QSH insulators 
by Nie et al^, Sun et al^ and Ma et al^ using first-principles calculations. 

The paper is organized as follows. In section [^the effective tight-binding (TB) model of 


TMD haeckelites are constructed from symmetrical analysis. In section |HI[ the symmetry 
and topology of the lattice model and the topological phase transition have been discussed. 
In section |IV[ the vorticity, merging or annihilation of Dirac cones have been demonstrated 
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by k-p model. In section |V| we provide a brief snmmary. 


II. DERIVATION OF TIGHT BINDING MODEL 

Monolayer of TMD in its top view, snch as M 0 S 2 , can be thonght as a honeycomb lattice 
with Mo and S on the A and B sublattice site, respectively. However, recently a kind 
of ordered square-octagonal distortion have been introduced in honeycomb TMD to form 
haeckelite structurd^^ES] graphene. Such TMD haeckelites have pairs of square- 

octagon composed of M and X atoms, noted as MX2-4-8 to be distinguished from usual MX 2 
with honeycomb lattice. In our previous workP^, the hrst-principles calculation for WS2-4-8 
indicates the four bands around the Fermi level are mainly contributed by orbitals of 
M atoms. The p orbitals from X atoms are quite far from Fermi level and can be safely 
ignored. In this point of view, the square-octagonal lattice can be further simplihed to be a 
rectangular lattice with only four M atoms per unit cell as shown in Fig. 


A. Hamiltonian without SOC Term 

We take orbital on each M atom as local orbital basis and note it as Wannier orbital 
Wn{r — Ri). Here n=1...4 indicates four nonequivalent M atoms in one unit cell as shown in 
Fig. and Ri is the lattice index. Such rectangular lattice has space group Pbam and the 
following symmetries are satished. 


V: inversion symmetry operator. Mi (M2) and M3 (M4) are related by V. 

rriz- mirror operator with mirror plane perpendicular to normal of the 2D lattice plane. All 
M atoms keep identity. 

{^*2x1^10}: This is a non-symmorphic operation. After C 2 rotation around x-axis (along a 
lattice), a glide of (|, |, 0) (in unit of lattice vectors) is applied. This operation relates 
Ml (M3) with M2 (M4). 
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Therefore, there are only three independent hopping terms in tight-binding (TB) approx¬ 
imation and they are: 


{tci|Trok2) 

(1) 

(tCi|i/ok3t) 

(2) 

{wx\Hq\w^ 

(3) 


Here Hq is the Hamiltonian of system without considering of spin-orbit coupling (SOC). 
The hopping parameters ti, ^2 and ts are marked in Fig. [T] also. All these parame¬ 
ters must be real because and Hq are real. After fourier transformation Honni{k) = 
^^(tc„(0)|Tro|'n'm(-R))e*^''^, we get k dependent Hamiltonian HQ^k) satisfying periodic con¬ 
dition Ho{k)=Ho{k -h G) with G being integer times of reciprocal lattice. A k-dependent 
gauge transformation HQ^k) —)■ U{k)HQ{k)U^{k) is applied to get a concise but not periodic 
form of HQ(k) in Eq. Q. This form is more convenient for deriving the effective k-p models 
used in section HVl 


Ho{k) 


^ 0 

7i 

72 

73 

ll 

0 

73 

it 

l2* 

73 

0 

ll 

CO * 

72^* 

7i* 

0 


7 l = 2ti cos ^ 

73 = 2^3 cos 


U{k) 


/ 1 0 0 0 \ 

0 0 0 

0 0 0 

\0 0 0 j 


(4) 


(5) 


4 






B. Spin-Orbit Coupling Terms 


To include SOC term, the degree of freedom of spin is introduced explicitly to expand 
the dimension of the Hamiltonian from four to eight. The spinor basis set is dehned as: 

|0i) = Wi(r) 0 i = (6) 

\(j)i+4) = Wi{r) ^ Xi * = (7) 

, where X'|-=(1,0)^ and X4,=(0,1)^ with T representing transpose. Correspondingly, all the 
above symmetrical operators should be adjusted: 

V 


ruz —)■ rUzC 


1 1 


22 


11 


C'o.tI ——0 ^ < Cos.! ——0 >62 


22 


and TR operator is dehned as T = —iayK with K is complex conjugation operator. Con¬ 
sidering these symmetrical constraints, we will derive the term including SOC Hgo- 

Mirror symmetry mz ■ Similar as the case in Kane-Mele modeP^ on honeycomb lattice 
without Rashba term, Hgo for present TMD haeckelite is also block diagonal in spin space 
because of symmetry. This can be proved as: 


Vi, j = 1- ■ - 4 

(0i+4|i7so|0j) = {mze~i'^‘"^(f)i+i\Hso\mze~'^'"'^^(t)j) 

(02+41 ) 0 

The derivation from the second line to the third line is due to = —icx^. 

Time-reversal symmetry T. and are related by T: 

(0i+4|i7so|0i+4) = {T(t)i\Hso\T4)j) 

= (00 77*010i) (8) 


i.e., Hll = . 

Inversion symmetry V and combined symmetry ■ T■ These two opera¬ 

tions do not hip spin. After an intuitive graphic analysis as illustrated in Fig. [T| one can 
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easily find that there are only two independent non-zero parameters: 


( 011 I ^ 2 ) ( 031 104 ) '^'^1 

(0l|-f^so|04) = (03|-^so|02) = ^^3 


( 9 ) 

( 10 ) 


Note that these two SOC terms must be pure imaginary numbers because 0i are real and 
Hgo contains an imaginary unit. We also note that (0i|h?^so|030 must be zero because V 
transforms it to its complex conjugation. 

Therefore, Hgo is always diagonal in spin space and VT guarantees spin degeneracy, i.e., 
Kramer degeneracy. We only need to deal with for the rest of this article. The concise 
matrix form of Hgo in Eq. (11) is achieved by the same gauge transform (Eq. ([^) performed 
to one spin channel of Hso{k). 


HVoik) = 


( 11 ) 


f 0 0 53 \ 

S* 0 s; 0 

0 0 (5i 

V<53* 0 5,* 0 / 

= 2iAi cos 
^3 = 2 iA 3 cos 

Now we get the total effective TB model H{k)=HQ{k) + Hso{k), but keep in mind that 
the periodic condition is sacrihced to get the conciseness in Eq. Q and Eq. ( |IT| , and 
a gauge transformation (Eq. ([^) to periodic form is needed when calculating properties 
concerning global phase such as Berry phase and related properties. 


III. SYMMETRY AND TOPOLOGY OF THE BANDS 


From the analysis of symmetry in the above section, we obtain a TB model as shown in 
Eq. (|^ and Eq. (11). We now focus on the band structure and the topology with different 
choice of related parameters. In real materials, these parameters are determined by details 
in crystal structure and choice of M and X elements in MX2-4-8 haeckelites.^ It should be 
noticed that only the sign(tif 2 t 3 ) makes sense. This can be easily understood if we can use 
some unitary transformation UHU"’" to change tj’s signs in pairs. For example, U=l 2 <8) 
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(a) H„ 


(b)H 

^ SO 


FIG. 1: 2D rectangular lattice with four sties per cell, (a) Hopping parameters ti, t 2 and G 
in Hq. (b) Finite parameters Ai and A 3 in Hgo with each arrow points from bra to ket. 


transforms H, G to —G, —G and U=Tz ® I 2 transforms ^ 2 , G to —h-, —G. The change in the 
sign(tit 2 t 3 ), whatever which sign(tj) changes, is equivalent to Ho{k) —)■ —Ho{k). Without 
loss of generality, in the rest of this paper we deal with 0 < fa < fi. For the case |fi| < Ifsl, 
the system can be thought as the case |fi| > with the interchange of x and y axes, which 
is apparent in Fig. For clarity, we list all possible band structure with 0 < fa < fi in Fig. 

The classification and topology of them will be discussed in the following. 


A. Symmetry of The Eigenstates 

Ho{k) at time-reversal invariant momenta (TRIM)l^can be solved analytically. The eigen 
energies, degeneracy and eigenvalue of symmetrical operators of the eigen states are listed in 
Table Little group at F of Pbam has only one-dimensional irreducible representations,^ 
and all symmetrical operations can be diagonalized within the space spanned by eigen wave- 
functions. 

At TRIM X( 7 r, 0, 0) or Y (0, vr, 0), Pbam has only 2D irreducible representations. Within 
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TABLE I: Eigen energies, their degeneracy and eigenvalues of symmetrical operations of 
Hoik) at four time-reversal invariant momenta (TRIM). 



Energy 

Folds 

V 


{C2y\Ty} 

r 

Ea = ^2 — 2(ti -I- is) 

1 

1 

-1 

-1 


= —t2 — 2(ti — ts) 

1 

-1 

-1 

1 


Eq = —t2 -I- 2(ti — ts) 

1 

-1 

1 

-1 


E^ = t2 + 2(ti ts) 

1 

1 

1 

1 

A 


2 

±1 

±1 

±1 


\/^2 + 4^3 

2 

±1 

±1 

±1 

Y 


2 

±1 

±1 

±1 


\/tl + Atf 

2 

±1 

±1 

±1 

s 

—t2 

2 

1,1 

±i 

±i 


t2 

2 

-1,-1 

±i 

±i 


the eigen states of HQ{k), they can be written as following: 


Dx (P) = -Dy(P) = 


Dx[C^A\\o \ = 


Dy{C2,\\\o\ = 



( 12 ) 

( 13 ) 

( 14 ) 


in which r* denotes the i-th pauli matrix, representing the space spanned by two-fold de¬ 
generate eigenstates of Ho{X) or Ho{Y). We can character symmetry properties of each 
eigen space by eigenvalues of symmetrical operators. For example, the two-fold degener¬ 
ate eigenstates at either X (or Y) have even and odd parities respectively since inversion 
symmetry V is represented as (or T 2 ), thus we note this eigen space with ±1 in Table 
Solutions at S'(7r, 7r,0) are special. Little group of Pbam at S has only one-dimensional 
irreducible representations. However, for {C'2x|^^0} and {<^2^1^10}, the eigenvalues of the 
two eigenstates with the same parity are ±L Thus, the TR operator connects them and 
makes them two-fold degenerate. 

























FIG. 2: Band structure of Ho{k) for all possible parameter choices with 0 < ts < ti. (a)-(f) 
and (g)-(j) show bands without (fs 7 ^ G) and with (^3 = ti) symmetry, respectively. 
Upper and lower panels are for ^2 > 0 and ^2 < 0 cases, respectively. Color of each band 
represents its representation under the conserved non-symmorphic rotational symmetry. 

Along F-X (F-Y), black (green) and red (yellow) represent +1 and -1 eigenvalue of 
{C 2 x\\\^} respectively. Including SOC term will drive the system into QSH 

insulator or normal insulator (NI) as indicated for different parameter choice. 


From Table |T| we can see that the order of eigen energies at X and Y do not depend on 
the choice of f*. Further, since the parities of the double degenerate states are always in 
pair, the band inversion at X and Y will not change the topology of the occupied bands. At 
S', the situation is similar, though the energy order depends on sign of ^ 2 - Therefore, only 
the band inversion at F can result in different topology of occupied bands. This is similar 
to large-band gap 2D topological insulator ZrTes.^^ 

The order of energies at F depends on explicitly as illustrated in Fig. Along F-X(Y), 
{C 2 x\l^ 0 } {{C 2 y\\\^}) symmetry is conserved and the eigenstates with the same eigenvalue 
of {C 2 x\\\^} {{C 2 y\\\^}) at different k points form a continuous band, which is identihed 
by different color in Fig. As long as Id < 2ti, the {C 2 x\W^} eigenvalues of the two lower 
bands at F are -1,-1, while those at X are 1,-1. There must be a level crossing along 
F-X path. However, if Id > 2fi , the eigenvalues of the two lower bands at F and X are 
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the same and so there is no crossing. But we’d like to remind that even without the above 
C 2 symmetry, as long as TR and V is conserved, such spineless 2D lattice model will have 
stable band crossing once band inversion happens.®^ C 2 symmetry constraints the crossing 
point on the path along F-X or F-Y, which is similar to the spinless graphene model.^ A 
detailed discussion about these crossing points will be carried out in section [TV| 

B. Topological Phase 

We now investigate whether a QSH insulating state can occur or not after SOC term is 
considered, which is similar to what happens in Kane-Mele model.^^^ The Z 2 number of 
system with both TR and inversion V symmetries is determined by the product of parities 
of occupied states at four TRIM.I^ SOC term opens band gap at the band crossing points 
as shown in Fig. [^a). As long as 1 ^ 2 ! < 2 max(fi,f 3 ), Z 2 number is 1 as shown in Fig. 

To verify this we calculate the evolution of hybrid Wannier function centers for occupied 
bands using wilson loop method.®^ The plot is shown in Fig. [^b). This illustrates that 
the winding number of two occupied bands is 1, indicating that the Chern number of spin 
up bands is 1. Due to TR symmetry, the spin down bands should have Chern number — 1. 
Therefore, the model with |t 2 | < 2 max(fi,f 3 ) and non-zero SOC parameter Aj is a QSH 
insulator state. For the case |f 2 | > 2 max(fi,f 3 ), there will be no band inversion and the 
product of parities is 1, i.e., the system becomes a trivial normal insulator (NI). 

Explicitly, we investigate the topological phase transition (TFT) between TI and NI 
driven by tuning one parameter t 2 - Such TFT from QSH insulator to NI in centrosymmetric 
system will experience a gapless state with a Dirac node at TRIM.^^^^^As illustrated in Fig.|^ 
the TFT can be thought as the bands inversion of A and C at F in case 0 > 0 or R and D 
in case ^2 < 0. Hence, at TFT point the eigen energies of Ho(T) + HgoiX) satisfy Ea = Ec 
or Eb = Eb with 

Ea = t2 — + ^ 3 )^ -|- (Ai -|- A 3 )^ (15) 

Eb = -t2 - 2y/(H-t3)2 + (Ai-A3)2 (16) 

Ec = -t2 + 2^{ti - hf + (Ai - A3)2 (17) 

Ed = 0 + 2 \/(fi -|- + (Ai -|- A 3 )^. (18) 

All of the corresponding eigen states are spin degenerate. It’s easy to hnd the critical value 
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FIG. 3: Bands and wannier centers’ evolntion with SOC. (a) shows bands in which black 
lines denote non-SOC bands and red lines denote SOC bands, (b) shows wannier centers’ 
evolntion of spin np electrons. The dashed red lines indicate reference lines. The 
parameters are set as ti = 1, G = 1-3, G = 0 . 8 , Ai = —A 3 = 0.04. 


of G: 

Oc = ± V (G + + (Ai + Aa)^ 

± -n/(G - hy + (Ai - As)^ (19) 

Since t 2 c depends on U and A* {i=l and 3), band inversion is inflnenced by both hopping 
effect and SOC in real materials. 

IV. DIRAC POINTS AND QUADRATIC BAND TOUCHING 

From the above discnssion, onr model can be looked as an analog of graphene on a 
rectangnlar lattice: a signihcant resemblance is that QSH state occnrs when SOC term 
opens gap at Dirac points. Since the lattice constants a and b are qnite close to each other 
in previonsly stndied TMD haeckelites,!^ we’d like to go fnrther to see what will happen 
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FIG. 4: Berry curvature and pseudospin distribution in momentum space around around 
band crossing or touching points, (a) shows the Dirac cones with opposite vorticity. (b) 
and (c) show the Dirac cones having the same vorticity when 6*4 symmetry is slightly 

broken and conserved , respectively. 

in a square lattice with (^4 symmetry, which is the case studied in Ref. [2l] and |23l We 
hnd that in square lattice the non-SOC band structure has a quadratic band touching at 
r, being different from the linear Dirac cones in graphene and rectangular lattice. As all 
interesting physics are related to the linear Dirac cone or the quadratic non-Dirac cone band 
touching, in this section we use effective k-p model to characterize them. Alternatively, 
instead of calculating the evolution of centers of hybrid Wannier functions,®^ we calculate 
the integral of berry curvature analytically by treating Hso as an inhnity small perturbation 
term to identify the topological invariant. 


A. Linear Dirac Nodes in Rectangular lattice 

As the non-trivial topology comes from band inversion of A, C or B, D bands around 
r, it’s natural to choose the states of A, C or B, D as the bases to establish k-p model to 
capture the essence of the underlying physics. We take the case of A, C band inversion in the 
following discussion and the B and D inversion case can be easily inferred from it. We expand 
Eq. (|^ and Eq. ( pA| to second order of k and then use the downfolding techniqud^^® to 
reduce the dimensionality in order to get the effective Hamiltonian in the bases of eigenstates 
I A) and |G). To keep the k-p model similar to the Kane-Mele model, we transform |A) and 
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\C) to the following two bases: 


( 20 ) 


M\ _ -kl) + |W4) 1^^ _ \W2) - |W3) 

' ^ ^ s/2 

Within the above basis set, the space inversion and TR operators are 

DAci'P) = R and Dac{K) = K (21) 

, respectively, both of which are the same as those in Kane-Mele model. 

After some approximation (neglecting second order of Aj), we get the k-p model in basis 
set of | 1 ) and | 2 ): 


HAc{k) = EAc{k) + 

F = 


f 2 ( 2 - f 2 - ) + Fk'^ 




•'2 


8{Fa 


Eb) 

A_ 


8{Fc — Fb) 
A_ 


+ 


+ 


A, 


+ 


A_ 


'^1 “ -^kyT2 + Jk^n (22) 


(23) 

(24) 


_Fa — Fb Fc — Fb Fa — Fb Fc — Fb_ 
where FAc{k) is an even function of k and A± is dehned as A± = 2(Ai ± A 3 ). 


If SOC is not present, i.e., A±=0, there will be two Dirac cones from Eq. 22 and the 
location of the two nodes are determined by solution of HAc{k) = 0 . k^ is found to 
be: 


= Cl n 


2ti — ^2 


(25) 


F + (2ti — ^ 2)78 

Si is the unit vector along E-X path. In the case ti < t^, k^ is along F-Y. Expanding the 


hamiltonian Eq. 22 around k^ or —k^ to the lowest order of A± and q (here a general k is 
±k^+q), we get a linear dispersion hamiltonian, i.e. a Dirac hamiltonian: 


H{±k'^ + q) = ±ak^qiTi - ^-^q2T2 ± Jk^T^^ = d ■ t 


a = 


2 fi — ^2 


+ 2F 


(26) 

(27) 


Here d={di, d 2 , ds) and r=(ri, T 2 , rs). Representation of VK is tiK and {VK)~^H{k)VK = 
H[k) guarantees the mass term (the term proportional to rs ) to be zero. Therefore, d 
becomes a 2-component vector. The two Dirac cones possess opposite vorticity dehned as 
sign(det(c}d/9g)). In this sense, any perturbation holds VK symmetry can only move the 
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positions of Dirac cones bnt can not annihilate them. This can also be easily seen from the 
textnre of psendospin, which is dehned by psendospin operator r, As shown in Fig. |^a), the 
psendospin of occnpied state rotating 27r along an anticlockwise loop snrronnding a Dirac 
cone gives a positive vorticity, while —27r gives a negative vorticity. Psendospin rotates 0 
along a loop containing the two Dirac cones and a trivial insnlating state will be obtained 
when these two Dirac cones meet together. 

In the presence of SOC term, the time reversal operator becomes —iayK which is different 
from K in spinless case. Hence complex conjngation symmetry is broken and a mass term is 
introdnced which drives the system to QSH state as discnssed above. Using the formnlaP^ 

= id ■ {dk,d X dkyd) (28) 


where d = d/\d\, we calculate the berry curvature around Dirac nodes: 


d^Ac{q) = 


at2Jk{ 


a2 


fql + tlqijA + nf 


( 29 ) 


Integral of Eq. (29) around one Dirac point gives ^sign(af 2 </), thus the Chern number in 


spin up channel is 1 (or —1). Due to TR symmetry, spin Chern number would be 1 (or —1). 


B. Quadratic Band Touching in Square Lattice with 6*4 Symmetry 

In this subsection, we will consider the square lattice with 6*4 symmetry. This leads to 
fs = ti and A 3 = —Ai. Firstly, we consider Ffo(fc) without SOC. As shown in Fig. and 
[^, there is gapless band touching at F and it is shown to be quadratic in the following 
discussion. This is due to the degeneracy of Eb and Ec at F and band inversion of Ea and 
Ec- If band inversion is absent, it becomes a trivial normal insulator as shown in Fig. 
and [^. These two cases are similar to zinc-blend HgTe with band inversion and CdTe 
without band inversion 

Although the non-trivial topology comes from the A, C band inversion, the k-p model 
based on bands A and C can not reproduce the band touching at Fermi level because it 
comes from the degeneracy of B and C bands. Therefore, we construct a k-p model in bases 
of \B) and |C). To keep the similar form as Haci we also transform \B) and \C) to two new 
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bases: 

_ |wi) -i\w2) - IW 3 ) +i\wi) 

^ ~ 2 

|2^ ^ A^i) - \w2)-i\w3) + IW 4 ) 

Thus, the space inversion and complex conjugation operators are written as: 


( 30 ) 

( 31 ) 


DbcCP) = -1 Dbc{K) = -inK 


(32) 


which is different from Eq. and is responsible for the quadratic band touching. We get 
the Hamiltonian defined in basis set consisting of |1) and |2): 


Hscik) — —t 2 



+ 



n + 


h 


kxkyT 2 + A.rg 


( 33 ) 


Similar to the last subsection, by expanding the Hamiltonian around band touching points 
to the lowest order of A± and q, we get two Dirac cones: 


HBci±k^ + q) 


2 


qm ± 



q 2 T 2 + A.rg 


= Cl 



^3 

h 


( 34 ) 

( 35 ) 


If C 4 symmetry is slightly broken, the two linear Dirac cones are slightly separated by 2 |fc^|. 
They are related by space inversion. However, because V = —1 the two Dirac cones have 
the same vorticity. The pseudospin texture in Fig. |^b) indicates the two Dirac cones have 
the same vorticity since pseudospin rotates dvr along an anticlockwise loop surrounding the 
two Dirac cones. This is the reason of quadratic band touching at T when 6*4 symmetry 
brings these two Dirac cones together. Such quadratic band touching can also be seen from 


Eq. 33 by taking tg = ti and Ag = — Ai required by C 4 , symmetry: 


Hscik) = Escik) + - kl)Ti + ^-^kxkyT2 + A_rg 


(36) 


where Escik) is an even function of k. 

Similarly, SOC breaks complex conjugation symmetry and introduces a non-zero mass 
and removes the above quadratic band touching degeneracy at T. This leads the system into 
insulator. To identify the band topology of this insulating state, the momentum distribution 
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of berry curvature around F can be calculated by Eq. (28): 


7 = 




4W16A2 + 


kr 


t\ co7{26) + ^ sin^(26') 


9 = arccos ( -7 

\k\ 


(37) 


It is shown in Fig. 4(c). The four-fold symmetrical distribution can be easily seen and 
understood from the 26 terms. The integral of this berry curvature gives 1 (or —1). This 
indicates the Chern number in spin up channel is 1 (or —1) and the time-reversal symmetry 
ensures this insulating state to be a QSH insulator. As shown in Fig. |^c), the two Dirac 
cones are now moved together to T to form a double Dirac cone,l^ and the pseudospin 
rotates dvr after adiabatic loop evolution surrounding it,^ which indicates the quadratic 
band touching is non-trivial after Hso is included. 


V. CONCLUSIONS 

In this paper, we derived an low energy TB model for 2D TMD haeckelites MX2-4-8 to 
uncover the physical mechanism controlling the topological phase transition. This square- 
like or rectangle lattice model contains only one orbital per site and four sites per unit cell. 
The band inversion is tuned by hopping parameters, which brings linear Dirac cones with 
the same or opposite vorticity. The annihilation or merging of them has been demonstrated 
through effective k-p model analysis and elaborated by pseudospin texture. SOC plays 
a critical role in opening band gap at the band crossing or band touching points, which 
drives the system into insulator. The hybrid Wannier function center evolution calculated 
by Wilson loop method, as well as the direct integral of Berry curvature, is used to identify 
the topology of insulating state. The QSH insulator can be obtained within reasonable 
parameter region and can be used to understand the band structures of TMD haeckelites 
MX2-4-8 with M=Mo, W and X=S, Se, Te. Our model has extended the QSH lattice model 
family to square or rectangular lattice. 
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